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ABSTRACT 
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Motivated by the substantial sensitivity of eddies in two-layer quasi- 
geostrophic (QG) turbulence models to the strength of bottom drag, this study 
explores the sensitivity of eddies in more realistic ocean general circulation 
model (OGCM) simulations to bottom drag strength. The OGCM results are 
interpreted using previous results from horizontally homogeneous, two-layer, 
flat-bottom, f-plane, doubly periodic QG turbulence simulations and new re- 
sults from two-layer B-plane QG turbulence simulations run in a basin ge- 
ometry with both flat and rough bottoms. Baroclinicity in all of the simula- 
tions varies greatly with drag strength, with weak drag corresponding to more 
barotropic flow and strong drag corresponding to more baroclinic flow. The 
sensitivity of the baroclinicity in the QG basin simulations to bottom drag is 
considerably reduced, however, when rough topography is used in lieu of a flat 
bottom. Rough topography reduces the sensitivity of the eddy kinetic energy 
amplitude and horizontal length scales in the QG basin simulations to bottom 
drag to an even greater degree. The OGCM simulation behavior is qualita- 
tively similar to that in the QG rough bottom basin simulations in that baro- 
clinicity is more sensitive to bottom drag strength than are eddy amplitudes or 
horizontal length scales. Rough topography therefore appears to mediate the 
sensitivity of eddies in models to the strength of bottom drag. The sensitiv- 
ity of eddies to parameterized topographic internal lee wave drag, which has 
recently been introduced into some OGCMs, is also briefly discussed. Wave 
drag acts like a strong bottom drag in that it increases the baroclinicity of the 


flow, without strongly affecting eddy horizontal length scales. 
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1. Introduction 


This study focuses on the impact of frictional bottom boundary layer drag (“bottom drag” here- 
after) on the statistics of the midocean eddy field, where eddies are defined as deviations from a 
time-mean. The focus on bottom drag is motivated by the substantial sensitivity of eddy statistics 
to bottom drag strength documented in numerous studies of flat-bottom quasi-geostrophic (QG) 
turbulence (Salmon, 1978, 1980; Haidvogel and Held, 1980; Larichev and Held, 1995; Ozgoékmen 
and Chassignet, 1998; Riviere et al., 2004; Arbic and Flierl, 2004; Thompson and Young, 2006; 
Arbic et al., 2007; Arbic and Scott, 2008; Straub and Nadiga, 2014). A consistent finding in such 
studies is that weak bottom drag leads to a vigorous inverse cascade yielding a strongly barotropic 
and energetic eddy field characterized by horizontal length scales significantly larger than the first 
baroclinic mode deformation radius (Lg). Computations of spectral kinetic energy fluxes made 
from satellite altimetry, idealized models, and realistic ocean general circulation model (OGCM) 
simulations (e.g., Scott and Wang, 2005; Scott and Arbic, 2007; Schlosser and Eden, 2007; Qiu et 
al., 2008; Tulloch et al., 2011; Arbic et al., 2013, 2014; Straub and Nadiga, 2014) suggest that an 
inverse cascade to larger spatial scales is ubiquitous in the surface ocean. Indications are, however, 
that the inverse cascade proceeds over a relatively narrow range of oceanic length scales. Accord- 
ingly, observations demonstrate that the oceanic mesoscale eddy field lies far from the weak-drag 
limit of flat-bottom QG turbulence. For example, Wunsch (1997) finds that oceanic eddies are not 
strongly barotropic — instead, the kinetic energy levels in barotropic and first baroclinic modes are 
comparable. Stammer (1997) finds that length scales of ocean eddies are not much greater than Lg 
— instead, they are only slightly greater. Arbic and Flierl (2004) and Arbic and Scott (2008) argued 
that the “moderate” drag regime of QG turbulence (in between the weak drag and very strong drag 


limits) compared best to observations. However, it must be noted that most of the geostrophic tur- 
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bulence studies above are highly idealized, typically assuming not only QG dynamics, but in some 
cases also assuming horizontal homogeneity, zonal mean flows, a flat bottom, f-plane dynamics, 
and/or a severe truncation of vertical resolution to two layers. In some studies, the stratification 
is further simplified to consist of two layers of equal depths, thus precluding examination of the 
effects of surface-intensified stratification. The question therefore arises as to whether the sensi- 
tivities to bottom drag seen in the simple QG models used in many previous studies would also 
arise in more complex models such as high-resolution ocean general circulation models. 

More realistic OGCMs have rough topography, non-zonal mean flows, the planetary B-effect, 
surface-intensified stratification, ageostrophic dynamics, many layers in the vertical direction (not 
just two), and stratification and mean flows that vary in the horizontal direction. Any one of these 
factors could alter the sensitivity of eddy statistics to bottom drag. For example, Briiggemann and 
Eden (2015) have demonstrated that the routes to energy dissipation associated with ageostrophic 
and quasi-geostrophic flows are qualitatively different, with the energy flux towards smaller scales 
in (O(1) Rossby number) ageostrophic dynamics and towards larger scales in geostrophic turbu- 
lence. Increased vertical resolution implies that a lesser fraction of the water column will directly 
feel the effects of bottom drag, such that the sensitivity of eddy statistics to bottom drag is likely 
to be impacted. Horizontal inhomogeneities in more realistic models provide a more realistic 
environment for eddy evolution, and this may also affect eddy statistics (Merryfield 1998). Ve- 
naille et al. (2011) examined horizontally homogeneous QG turbulence simulations with a surface- 
intensified stratification, several layers in the vertical direction, imposed mean flows that project 
onto higher vertical modes, non-zonal mean flows, and the planetary B-effect. Similar to earlier 
studies, which often did not include many of these effects, they also found a strong sensitivity of 
the model eddy field to bottom drag strength (see their Table 2). Topographic effects, however, 


were not considered in their study, whereas it is well know that topography can profoundly influ- 
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ence the eddy field (Rhines 1970, 1977; Treguier and Hua 1988; Treguier and McWilliams 1990; 
Dewar 1998; Sinha and Richards 1999; LaCasce and Brink 2000; Benilov et al. 2004; Hurlburt et 
al. 2008; Thompson 2010; Boland et al. 2012; Venaille 2012; Chen and Kamenkovich 2013; Aber- 
nathey and Cessi 2014; Stewart et al. 2014; Chen et al. 2015). Some of these topographic effects 
involve small vertical length scales and are thus poorly represented in ocean general circulation 
models, which typically concentrate vertical resolution near the surface. One result of particular 
interest from the studies mentioned above is that topography can facilitate a downward transfer of 
energy (Venaille 2012). Note that at (forced-dissipated) statistical equilibrium, this need not imply 
a strong bottom intensification of kinetic energy because kinetic energy is continually input by the 
forcing and abyssal energy is removed by bottom friction. 

Two additional factors typically absent in idealized studies, but that might also influence ocean 
eddy statistics, are internal lee waves and topographic blocking (together referred to as “wave 
drag” hereafter). Interest in wave drag, as a contributor to the oceanic energy budget and a poten- 
tially important addition to ocean model dynamics, has grown rapidly in recent years. The internal 
lee wave contribution to wave drag is the momentum flux due to wave generation over certain 
topographic length scales. The topographic blocking contribution to wave drag occurs when the 
streamline is parallel to the seafloor, and characterizes the hydraulic effects, low-level breaking, 
vortex shedding, flow separation, and low-level jets (Baines, 1995) that occur when flow impinges 
upon a topographic feature. Using a closure first developed by Garner (2005), Trossman et al. 
(2015) compared predictions of dissipation profiles in the Southern Ocean with microstructure 
profiler observations, and argued that the topographic blocking contribution to wave drag domi- 
nates the dissipation in the bottom 1000 meters. Trossman et al. (2013, 2016) found more than 0.4 
TW of low-frequency mechanical energy dissipation associated with the combination of internal 


lee wave generation/breaking and topographic blocking in a model run with the Garner (2005) 
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wave drag parameterization. Nikurashin and Ferrari (2011), Scott et al. (2011), and Wright et al. 
(2014) all estimate that breaking internal lee waves dissipate at least 0.2 TW of low-frequency 
mechanical energy, comparable to the amount (0.1 — 0.2 TW) of dissipation estimated to occur 
via bottom drag (Sen et al. 2008; Arbic et al. 2009; Trossman et al. 2013; Wright et al. 2013; 
Trossman et al. 2016). Internal lee waves have also been found to be important in the momentum 
and vorticity budgets (Naveira-Garabato et al. 2013). 

Wave drag parameterizes ageostrophic effects and can be thought of as distinct from form drag. 
The latter is a correlation between bottom pressure and topographic slope. It can be thought of 
in terms of geostrophic dynamics and is known to be particularly important in the Antarctic Cir- 
cumpolar Current (ACC). Without form drag, closing the zonal momentum budget of the ACC in- 
volves either very large bottom drag or very large circumpolar transports (e.g., Olbers et al. (2004) 
and many others). In this context, recent work has explored the combined roles of bottom drag 
and topography in ACC settings. Various studies (Hogg and Blundell 2006; Nadeau and Straub 
2012; Nadeau and Ferrari 2015) have shown circumpolar transport to increase with bottom drag. 
This can be easily understood in the strong drag limit of the quasi-geostrophic equations. In this 
limit, abyssal velocities are weak, so that the bottom layer streamfunction (equivalent to pressure 
in quasi-geostrophy) becomes near-constant. As such, form drag is diminished and circumpolar 
transport is increased. Primitive equation models also show transport to increase as bottom drag 
coefficients are made large, although it is likely that the degree to which circumpolar transport 
depends on the bottom drag may be related to the complexity of the bottom topography and may 
be less than implied by these idealized studies (e.g., Nadeau et al., 2013; Nadeau and Ferrari, 
2015). 

In this study, we compare eddy statistics across realistic high-resolution OGCM simulations 


with varying strengths of bottom drag. For simplicity, the OGCM simulations analyzed here do 
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not include tides. In order to tease out the sensitivities very clearly, we vary the bottom drag 
coefficient Cg by a large factor (~ 500). Estimates of Cg values in the ocean vary by much less 
than that. Observations of boundary layer turbulence suggest Cz values of about 0.0025, with an 
uncertainty of a factor of about 3 in either direction (Weatherly, 1975; Trowbridge et al., 1999; 
Trowbridge and Elgar, 2001; Feddersen et al., 2003). We focus here on the statistics that Arbic 
and Flierl (2004) and Arbic and Scott (2008) focused upon — eddy baroclinicity or “vertical struc- 
ture,” eddy horizontal length scales, and eddy amplitudes — in their examination of the impact of 
bottom drag on two-layer flat-bottom QG turbulence. We compare the OGCM sensitivities to bot- 
tom drag with the sensitivities seen in previous studies of horizontally homogeneous, two-layer, 
flat-bottom, f-plane, doubly periodic QG turbulence, and the sensitivities seen in new two-layer, 
B-plane QG basin turbulence runs with both flat-bottom and rough-bottom conditions. Compar- 
ison of the multi-layer OGCM versus two-layer QG simulations may potentially shed light on 
the importance of ageostrophic effects and vertical resolution. Comparison of the horizontally 
homogeneous versus basin QG simulations may shed light on the importance of flow inhomo- 
geneities. Comparison of the flat-bottom versus rough-bottom QG basin simulations illuminates 
the importance of rough bottom topography in setting the sensitivity of eddying flows to bottom 
drag strength. Motivated by the growing interest in wave drag, this paper will briefly discuss the 
impact of wave drag upon eddy statistics by examining the OGCM simulations run with wave drag 
in Trossman et al. (2013, 2016). We note that Hurlburt and Hogan (2008) also did simulations 
of an OGCM with varying values of bottom drag. They used an OGCM (the Naval Research 
Laboratory’s Layered Ocean Model, NLOM) that is in a realistic domain, albeit with a number 
of simplifications relative to HYCOM. Hurlburt and Hogan (2008) focused on the response of 
western boundary current dynamics to bottom drag rather than on the impact of bottom drag on 


the inverse cascade of geostrophic turbulence. 
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The present paper is organized as follows. We first describe the high-resolution OGCM sim- 
ulations, carried out in both Atlantic Ocean and global domains assuming different bottom drag 
parameter values, and in the global domain with and without wave drag. We then describe the 
B-plane QG basin simulations, carried out in a midlatitude double gyre setting —with and without 
rough topography and assuming different values for a bottom drag parameter. We also briefly dis- 
cuss the setups for the Arbic and Flierl (2004) and Arbic and Scott (2008) two-layer, flat-bottom, 
horizontally homogeneous QG turbulence simulations that we will use here. We next describe 
various diagnostics used to measure the baroclinicity, amplitudes, and horizontal length scales of 
midocean eddies. Finally, we discuss the impact of bottom drag on eddy statistics in the QG and 
OGCM simulations, and the impact of wave drag on eddies in OGCM simulations. The diagnos- 
tics and results sections use some current meter observations and satellite altimeter products for 
comparison to the OGCM results. We end with some concluding remarks about the implications 


of this study. 


2. Model configurations 


The nominally 1/12° and 1/25° HYbrid Coordinate Ocean Model (HYCOM; Bleck, 2002; 
Chassignet et al., 2003; Halliwell, 2004) simulations are on a tripole Mercator grid and have 
32 hybrid layers in the vertical direction. HYCOM smoothly transitions between different vertical 
coordinates, depending on the relative strengths of the coordinates in different oceanic regimes 
(Griffies et al. 2000; Chassignet et al. 2006). The vertical coordinates are isopycnal in the subsur- 
face open ocean, z-level in the open ocean mixed layer, and terrain-following in shallow regions. 
The performance of HYCOM without wave drag has been evaluated extensively in the North At- 
lantic (Xu et al., 2016, and several references therein), in the North Pacific (Kelly et al. 2007), in 


the Indian Ocean (Srinivasan et al. 2009), and across the entire World Ocean (Chassignet et al. 
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2009; Thoppil et al. 2011). The performance of HYCOM with wave drag has been evaluated by 
Trossman et al. (2016) across the entire World Ocean. 

We now discuss the vertical and horizontal eddy viscosity parameterizations in HYCOM. The 
K-Profile Parameterization (KPP; Large et al., 1994) yields relatively strong vertical mixing in 
the mixed layer, with a smooth transition to weaker vertical mixing below. Background mixing is 
typically used in deep water with an assumed Prandtl number of three so that the vertical viscosity 
is a factor of three larger than the vertical diffusivity. Shear instability mixing is typically used 
in the mixed layer with an assumed Prandtl number of one. The horizontal viscosity includes the 
maximum of a Smagorinsky (1993) parameterization or Laplacian term with an additional bihar- 
monic term (Chassignet and Garraffo 2001; Chassignet and Marshall 2008). Horizontal viscosity 
smooths out subgrid-scale noise. Here, “horizontal” means following a vertical coordinate layer. 

For the global 1 /25° runs, we begin with a simulation that is spun-up using 1.125° x 1.125° Eu- 
ropean Centre for Medium-Range Weather Forecasts (ECMWF) Re-Analysis (ERA-40) monthly 
mean forcing over 1978-2002 (Kallberg et al. 2004; Uppala et al. 2005), supplemented with higher 
frequencies. Six-hourly anomalies with respect to monthly means from the 2003 fields of the Navy 
Operational Global Atmospheric Prediction System (NOGAPS; Rosmond et al., 2002) are added 
to the ERA-40 climatological wind forcing. The six-hourly winds are used during every model 
year in this way. 

The global 1/25° HYCOM simulation described above is first spun-up from rest for thirteen 
years using a value of the bottom drag coefficient (Cy = 2.5 x 10°) that is designated as “mid” 
hereafter. The mid Cy value is the reference, or “control,” value used in most HYCOM simulations. 
For legacy reasons, there is an assumed background tidal velocity (see, e.g., Willebrand et al., 
2001) of 5 cms~! for the first one and one-half years. The background tidal velocity is reduced to 


2cms~! for the next two and one-half years and 0 cm s~! thereafter. Starting at the end of year 12, 
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this HYCOM simulation is further integrated in two different configurations. One configuration 
is run for an additional 5 years with Cz = 2.5 x 107! (designated “strong” hereafter). The other 
configuration is run for an additional 4 years with wave drag and the mid value of bottom drag 
(Trossman et al. 2013, 2016). Daily averages of vertical velocity profiles at select locations, daily 
averages of sea surface heights, and bi-monthly averages of all other diagnostic model output are 
saved during the final year (year 13 for the mid drag value, year 17 for the strong drag value, and 
year 16 for the wave drag simulation). Because all of the results in this paper are computed from 
years that are well beyond the years in which there is a legacy background tidal velocity, the legacy 
tidal velocity does not affect any of our conclusions here. 

Only the 1/12° Atlantic configuration is run with the weak value of the bottom drag coeffi- 
cient (C7 = 5.0 x 10~*). The main reason for this is that simulations with the weak bottom drag 
coefficient require a very small baroclinic time step, making a global weak drag simulation pro- 
hibitively expensive. The 1/12° Atlantic simulation is first spun-up from rest for twenty-three 
years with a mid bottom drag coefficient (Cz = 2.5 x 1073). Sixteen spin-up years have a 5 cm 
s~! background tidal velocity and another seven years have no background tidal velocity. This 
simulation is then integrated for an additional 4 years with the weak value of the bottom drag 
(Cy = 5.0 x 1074). Daily averages of vertical velocity profiles at select locations, daily averages 
of sea surface heights, and monthly averages of all other diagnostic model output are saved during 
the final year (year 23 for the mid drag simulation and year 27 for the weak drag simulation). 
Table | presents the Cg values as well as the barotropic and baroclinic time steps of the HYCOM 
simulations analyzed in this paper. Note that both the weak and strong bottom drag runs require 
much smaller baroclinic time steps than the mid strength bottom drag (or control) runs. The wave 


drag simulation also requires a smaller time step. 
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The flat bottom QG B-plane basin model configuration used here is taken directly from Straub 
and Nadiga (2014). It has a uniform horizontal grid with Ax ~ 7.8 km, or about four grid points 
per deformation radius, Lz, here taken to be 30 km. The number of grid points is 512 x 512. The 
upper and lower layer thicknesses are set to be 1000 and 3000 meters, respectively. A double gyre 
(i.e., sinusoidal) zonal wind-stress forcing is applied to the upper layer potential vorticity equation. 
Biharmonic friction is added to damp enstrophy. A version of free slip conditions appropriate for 
biharmonic dissipation is applied; specifically, both vorticity and its Laplacian are set to zero 
at the horizontal boundaries. A Rayleigh (linear Stommel bottom) drag is applied to the lower 
layer only. The QG basin simulations analyzed here differ only in their bottom drag coefficient 
(rog = 8.0x 10719 s~!, rag = 8.0 x 1078 s~!, and rag = 8.0 x 10~° s~! are used, with the middle 
value taken as the nominal value) and in their bottom boundary condition (flat bottom and rough 
bottom topography). The rough topography used is taken from the North Atlantic region of the 
Smith and Sandwell (1997) bathymetric product. We want a topography that is rough, but is not 
rough at the model grid scale, as the latter would lead to numerical noise. In order to achieve this, 
we perform a two-dimensional interpolation of the Smith and Sandwell (1997) topography to a 
uniform 128 x 128 grid in the region bounded by 18.0 — 54.1°W, 7.3 — 43.4°N, and then perform 
another interpolation to the model’s 512 x 512 grid within the same domain. The bathymetry used 
in our rough bottom QG simulations is shown in Fig. 1. The figure shows the Mid-Atlantic Ridge 
cutting through the domain from the upper-right towards the lower-left corners. Note that our 
topography violates the QG assumption that the bottom layer depth variations are much less than 
the total depth. We also note that, as in most QG double gyre simulations, the formal requirements 
that BL/fo and €/fo be small are also violated, for linear meridional gradient in the Coriolis 
parameter B, topographic horizontal length scale L, Coriolis parameter fo, and relative vorticity 


¢. The time-averaged total energies are saved for each of the six QG basin model configurations, 
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following an initial spin-up sufficient to allow for energy levels to equilibrate. Daily output is 
saved for the ensuing final 135 days beyond the initial spin-up. 

The horizontally homogeneous, two-layer, flat-bottom, f-plane, doubly periodic QG results are 
taken from Arbic and Flierl (2004) and Arbic and Scott (2008). A linear bottom drag was used 
in the former paper while a quadratic bottom drag was employed alongside a linear drag in the 
latter paper. Arbic and Scott (2008) demonstrated that the impacts of bottom drag strength on the 
vertical structure, amplitude, and horizontal length scales of eddy kinetic energy are qualitatively 
similar whether linear or quadratic bottom drag is used; however, the sensitivity to drag is reduced 
when the drag is quadratic. The horizontally homogeneous QG results are run in a doubly periodic 
domain, with an imposed, baroclinically unstable mean flow meant to mimic the flows in a mid- 
ocean gyre. Equilibration results when the energy extracted by eddies from the mean flow is 


balanced by the energy dissipated by bottom drag. 


3. Diagnostics and observations 


For the most part, we compare our various model simulations with each other. However, we will 
also compare the sea surface height (SSH) variance, the geostrophic surface kinetic energy (SKE), 
and the vertical structure of the kinetic energy (KE) of the OGCM simulations with observations. 
The “observed” geostrophic SKE and SSH variance are taken from satellite altimetry products. 
To make the observations comparable with our model output, a mean SSH product (Andersen 
and Knudsen 2009; Andersen 2010) is added to the SSH anomalies from satellite altimetry before 
computing the observed geostrophic SKE and SSH variance. Geostrophic SKE is computed from 
the SSHs using a nine-point stencil according to the method outlined in Arbic et al. (2012). The 


model’s SSH variance and geostrophic SKE are calculated from daily averaged model output. 
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KE profiles at current meter locations (taken from the Global Multi-Archive Current Meter 
Database)! will be compared to the output of our global HYCOM simulations. The current meter 
velocities were filtered using a Butterworth filter with half-power of 3 and a daily cutoff period 
to eliminate tides and other higher frequency motions that are not present in the daily-averaged 
model output. We show the average vertical profile of the KE computed over the locations where 
current meter observations of at least a month duration exist. We place the KE at each horizontal 
location into 500 meter depth bins in the upper 4000 meters because 500 meters is a typical vertical 
resolution of abyssal layers in HYCOM; the vertical spacing between current meters on a typical 
mooring is of the same order of magnitude. 

We measure the vertical structure, or baroclinicity, of eddy KE in two ways: as the ratio of the 
baroclinic to barotropic KE (KEgc to KEgr) and as the ratio of near-surface to near-bottom KE. 
Here, KEgr is the kinetic energy of the depth-averaged flow, and KEgc is the kinetic energy of the 


deviations from the depth-averaged flow. For QG, 


AW + AoW 

SO is 

wee H, +> (1) 
Wi— Wo 
= /H,Ao———_ 
WBC 1 2 + Hy? 


where H and yy are the layer thickness and streamfunction in the upper layer, and Hz and > are 
the layer thickness and streamfunction in the bottom layer. Arbic and Flierl (2004) found the KEgc 
to KEgr ratio to be a more useful diagnostic for quantifying baroclinicity in weak bottom drag QG 
turbulence simulations, while the surface-to-bottom KE ratio was more useful in the strong drag 
limit. Only in our Atlantic simulations do we quantify baroclinicity using KEgc/KEgr. In both 
our Atlantic and global simulations, we use the top 100 meters and bottom 500 meters to represent 


the near-surface and near-bottom ocean; we will refer to the ratio of the two as KE;op100 /KEbotso0- 


'See: http://stockage.univ-brest.fr/“scott/GMACMD/updates.html (Scott et al. 2010). These observations were quality controlled by Timko et 


al. (2013) for effects such as blow-over. 
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This choice is made because the surface mixed layer is typically on the order of 100 meters thick, 
while the bottom two layers together in HYCOM are typically about 500 meters thick. When 
calculating the ratios, we omit all grid points where the water is shallower than 500 meters. When 
tabulating the area-averaged KEgc/KEgr and KE;op100 /KEpo1so0 ratios, we also omit all grid 
points within 30 indices of the coasts because in such locations there can be infinitesimal layer 
thicknesses that lead to finite transports but very large values of KE. 

Eddy horizontal length scale diagnostics are also computed. As in the doubly periodic QG 
turbulence simulations of Arbic and Flierl (2004) and Arbic and Scott (2008), we examine the 
length scales Lxg of eddy SKE and Lz of eddy barotropic KE. The HYCOM eddy SKE length 
scales are computed assuming a geostrophic streamfunction, y = gn/f, where 7 is the daily- 
averaged SSH, g = 9.806 ms”? is the acceleration due to gravity, and f is the Coriolis parameter. 
The SKE length scale is 


ee = (2) 


J Exe(K)dk 


where x is the isotropic horizontal wavenumber and Exg(K) = k*|@|? is the geostrophic SKE 


Lee = 


spectrum, where * denotes a Fourier transform. The QG model’s eddy SKE length scales are 
computed from the upper layer’s streamfunction. The QG model’s eddy length scales associated 
with KEgr are calculated similarly, but using Egr = |dgr|* +|?gr|? in place of Exg. Because it 
suffices to show the HYCOM KEgzr fields for the conclusions we draw about Lar, HYCOM Ler 
are not calculated. The two-dimensional Fourier transforms above are calculated using data from 
20° x 20° regions. Using output from our HYCOM simulations, y is interpolated onto a uniformly 
spaced (= 7.8 km) latitude-longitude grid. The temporal mean and spatial trends within each box 
were removed for the HYCOM simulations. For the QG basin simulations, the temporal mean 
trend within each box was removed; no interpolation was necessary since these data were output 


on a uniformly spaced grid. 
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Because of their relevance to interpreting the differences between the simulations with varied 
bottom drag strengths and the simulations with wave drag, we describe the bottom drag and wave 
drag contributions to the KE equation. This KE equation can be written as in Trossman et al. 


(2013) 


PRE st + PRE adv = fF press + Pinput — © out put + CKE—PE (3) 


Here, Px is the time derivative of the globally integrated KE, Pxg ady is the KE change due to 
advective fluxes across the sea surface, Ppress is the divergence of KE associated with pressure 
differentials at the sea surface, Pini: 1s the wind energy input, Pourpur 18 the sum of all dissipative 
terms such as bottom drag and wave drag (see below), and Cxg-_s,pe is the conversion rate of 
KE to PE. Because of the form of the wave drag parameterization in the momentum equations 
(Trossman et al. 2013, 2016), it can be thought of as a linear bottom boundary layer drag with a 
spatially varying coefficient (rg;ag). The energy dissipation rate due to quadratic bottom boundary 
drag is given by Taylor (1919) 


Dgp = PoCa|up|°. (4) 


The energy dissipation due to a combination of topographic blocking and internal lee wave drag is 


given by Trossman et al. (2013) 
Dwp = PoFdrag|4a\”- (5) 


Here, Po = 1035 kg m7? is the average density of seawater with respect to 2000 dbar; Cy is 
the quadratic drag coefficient; u, is the velocity averaged over the bottom Hgp meters and uy is 
the velocity averaged over the bottom Hwp meters, with |-| indicating a magnitude; Ydrag 1S a 
positive-definite decay rate times a vertical length scale, computed from uy and a power spectrum 
associated with the underlying topography; Hgp = 10 meters is the height range above the seafloor 


(up to the surface if shallower than 10 meters) over which quadratic bottom drag is applied in the 
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model; and Hwp = 500 meters is the height range above the seafloor (up to the surface if shallower 


than 500 meters) over which wave drag is applied in the model. 


4. Results 


Using horizontal eddy length scales, KE budget terms, geostrophic SKE, SSH variances, and ra- 
tios of KEgc to KEgr and near-surface to near-bottom KE, we will evaluate the impact of bottom 
drag strength on HYCOM and QG f-plane basin dynamics. We compare sensitivities in our HY- 
COM and QG basin simulations with results based on simpler two-layer, flat bottom, horizontally 
homogeneous QG turbulence studies. We also compare the SSH variance, geostrophic SKE, and 
vertical structure of KE from our HYCOM simulations with observations to assess the degree to 
which the bottom drag strength (Cy) is important in maintaining realistic eddy statistics in these 
simulations. We finish this section by examining how eddy statistics are altered upon addition of 


wave drag, using the metrics described above. 


a. SSH variance and geostrophic SKE 


The area-averaged geostrophic SKE in the HYCOM simulations, which have realistic 
bathymetry, is relatively insensitive to bottom drag strength, being only slightly increased with 
larger bottom drag strength and slightly decreased with smaller bottom drag strength (Figs. 2a-d; 
Table 2). This contrasts with previous studies of two-layer flat-bottom doubly periodic QG turbu- 
lence simulations (e.g., Arbic and Flierl 2004; Arbic and Scott 2008) for which the sensitivity is 
much greater.” 

SSH variance shows a somewhat larger sensitivity (Fig. 3; Table 2). For example, the strong 


bottom drag simulation shows greater SSH variance in the Gulf Stream Extension than is the case 


The geostrophic SKE is larger in each of the HYCOM model simulations than in AVISO (Fig. 2e). This is due to a known deficiency of energy 
in the AVISO product (e.g., Chelton et al., 2011). 
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for the control run (Figs. 3a and 3c). This is also true in the intensified jet regions outside that 
of the Gulf Stream as well. Conversely, the weak bottom drag run shows less SSH variance in 
energetic currents (Figs. 3b and 3d). 

We infer that the changes in SSH variances shown in Fig. 3 are due to increased near-surface 
eddy-driven mixing in the strong bottom drag simulations. Radko et al. (2014) postulates that 
eddy-driven mixing increases with shear, and we find evidence that the near-surface shear increases 
with drag coefficient (see the discussion of Fig. 4 below). Furthermore, the ageostrophic flow is 
affected through the curl of the wind stress, mostly in regions with intensified jets (not shown). 
We surmise that there are alterations in baroclinic instability due to differences in an inferred 
conversion rate between kinetic and potential energy change with varied bottom drag strength. 
Trossman et al. (2013, 2016) argued, making reference to (3), that the conversion rate between 
kinetic and potential energy must change when wave drag is included and the same energetics 


argument holds for our experiment with increased bottom drag strength. 


b. Vertical structure of the kinetic energy 


The vertical structure of KE in our strongly damped HYCOM simulations is qualitatively con- 
sistent with that seen in idealized QG turbulence simulations, but agrees poorly with observations. 
Table 3 demonstrates that the ratio of KE in the upper to lower layers is greatly increased in the 
strong drag HYCOM experiment, as would be anticipated from strong drag horizontally homo- 
geneous QG turbulence results (Arbic and Flierl 2004; Arbic and Scott 2008). Fig. 4b shows 
KE profiles for the low-passed observations and the global 1/25° strong- and mid-strength bot- 
tom drag simulations. Data are temporally averaged at each location in the Global Multi-Archive 
Current Meter Database and then averaged over all locations shown in Fig. 4a. Strong bottom 


drag renders a more baroclinic, surface-intensified flow. The KE from the strong drag simulation 
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(red curve in Fig. 4b) is greatly reduced near the seafloor and less so at shallower depths. The 
poor comparison between the strong-drag run and observations suggests that the real ocean is not 
in a strong-drag regime, consistent with the conclusions of Arbic and Flierl (2004) and Arbic and 
Scott (2008). Baroclinicity in the weak- versus mid-drag HYCOM simulations also behaves in a 
qualitatively similar way to what is observed in horizontally homogeneous QG turbulence (Arbic 
and Flierl 2004; Arbic and Scott 2008). Table 3 suggests that the Atlantic weak-drag simulation 
is more barotropic (less surface-intensified) than the mid-drag simulation. 

We next consider geographical distributions of baroclinicity. Figure 5 shows maps of 
KE;op100 /KEbotsoo for the global mid- and strong-drag HYCOM simulations and Fig. 6 shows 
maps of the same quantity for the Atlantic mid- and weak-drag HYCOM simulations*). The loca- 
tions where KE shows strong baroclinicity in the global maps of Fig. 5 tend to be within 40° of the 
equator or confined within bands in the Southern Ocean. Fig. 5 indicates that the number of grid 
points that are highly baroclinic is greater in the strong-drag simulation than in the mid-drag sim- 
ulation, consistent with expectations from horizontally homogeneous QG turbulence simulations. 
In the weak drag simulations, baroclinicity is considerably reduced (compare Figs. 6a and 6b). 
Overall, baroclinicity of the KE in HYCOM behaves qualitatively as one might expect from from 
idealized, flat-bottom, horizontally homogeneous QG turbulence simulations: the flow becomes 
distinctly more barotropic with weak drag and more baroclinic with strong drag. An important 
difference from this classical picture is that surface and barotropic KE are individually less sensi- 
tive than is the case in classic studies of QG turbulence. This can be seen by inspection of Fig. 7, 


which displays KEgr in the North Atlantic for the global and Atlantic HYCOM simulations with 


3We did not save the total or baroclinic component of the KE for the 1/25° global simulations due to the large hard disk space requirements 
needed to save these fields. The combination of Figs. 4b and 5 with Table 3 are sufficient to demonstrate that the flow becomes more baroclinic 


with stronger bottom drag. 
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varying bottom drag strength. Although KEgr is weaker when bottom drag is stronger (Fig. 7), 
this dependence is much less pronounced than is the signal as seen in baroclinicity (Figs. 5 and 6). 

Our QG basin simulations allow us to examine the impacts of rough topography and lateral inho- 
mogeneities on eddy statistics in QG flow. Fig. 8 displays the baroclinicity (quantified with both 
of the measures discussed earlier), as well as the surface and barotropic eddy horizontal length 
scales, in the QG basin simulations (with both rough and flat bottom topography), the previously 
reported horizontally homogeneous, two-layer, flat-bottom, f-plane, doubly periodic QG simula- 
tions of Arbic and Flierl (2004) and Arbic and Scott (2008), and the OGCM simulations. The 
abscissa of Fig. 8 represents the nondimensional friction strength, as defined by Arbic and Scott 
(2008) for the doubly periodic simulations, and defined by the ratio of the friction value to the 
nominal, or “control” value, for the QG basin and OGCM simulations. The QG basin simulations 
show that increased bottom drag leads to a more baroclinic flow, as expected (see blue curves 
in Figs. 8a-b), and in qualitative consistency with the QG turbulence results shown in Figs. 8a-b 
(black curves). Also as expected, overall there is less KE in the QG basin simulations when bottom 
drag strength is increased (Table 4).4 However, the sensitivity of baroclinicity and eddy energy to 
bottom drag strength is greatly reduced from what is seen in the horizontally homogeneous QG 
turbulence results, especially when rough topography is introduced into the QG basin simulations 
(e.g., compare the squares-solid blue curve to the diamonds-dot-dashed blue curve relative to the 
black curves in Figs. 8a-b, and the much greater sensitivity in Table 4 for the flat versus rough 
bottom simulations). This reduced sensitivity relative to horizontally homogeneous QG turbu- 
lence results is also seen in the HYCOM simulations over areas of rough topography, e.g., over a 


sub-domain of the North Atlantic (between 59.3° — 39.3°W and 19.6° — 39.6°N) close to the one 


4The eddy kinetic energy is only at a level near that of observations when the bottom drag coefficient lies in a particular range, but this range is 


considerably broader when rough topography is present than when a flat bottom is employed. 
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shown in Fig. 1 (red curves in Figs. 8a-b). It seems clear that rough topography accounts for 
much of the discrepancy between our HYCOM simulations and expectations from classic studies 


of flat-bottom QG turbulence. 


c. Surface eddy horizontal length scales 


We next consider eddy horizontal length scales. In our HYCOM simulations, length scales Lxg 
associated with SKE are fairly insensitive to bottom drag strength (Fig. 8d; Table 5). Although 
we did not explicitly calculate a length scale for the KE in the barotropic mode of our HYCOM 
simulations, visual inspection of Fig. 7 suggests that it too is relatively insensitive to bottom drag 
strength. In contrast, the surface eddy horizontal length scales increase more dramatically with 
reducing drag strength in the weak-drag limit of the horizontally homogeneous, two-layer, flat- 
bottom, f-plane, doubly periodic QG turbulence results of Arbic and Flierl (2004) and Arbic and 
Scott (2008), as can be seen in Fig. 8d. The increase in surface length scales in these previous 
simulations is mainly due to an increase in the barotropic length scale (Fig. 8c). 

To investigate a possible reason for the weak sensitivity of HYCOM eddy horizontal length 
scales to bottom drag relative to flat-bottom, horizontally homogeneous QG turbulence results, we 
compare eddy length scales from our QG basin simulations with and without rough topography. 
We consider eddy length scales associated with barotropic KE (Fig. 8c) and surface, or upper 
layer, KE (Fig. 8d). As with the HYCOM simulations (red curves in Figs. 8d), there is no general 
trend for the eddy length scales as a function of bottom drag strength in our rough bottom QG basin 
simulations. However, the eddy length scales in the flat-bottom QG basin simulations behave more 
like the previous flat-bottom doubly periodic QG turbulence results — both barotropic and surface 
eddy length scales increase greatly as drag is weakened in the weak drag limit. Overall, our results 


suggest that rough topography reduces the sensitivity of eddy horizontal length scales to bottom 
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drag. This insensitivity can be visualized through examination of snapshots of the upper layer 
streamfunction, shown in Fig. 9, for the QG basin simulations. The flat bottom simulations show 
large qualitative differences as drag strength is altered. With rough topography, this sensitivity is 
markedly reduced. In addition, we note that the presence of topography matters less to the surface 
streamfunction when the drag is strong. For instance, the streamfunctions for the simulations 
with strong-drag in flat- and rough-bottom configurations (Figs. 9c and 9f, respectively) are more 
similar to each other than are the streamfunctions for the simulations with mid- or weak-drag in 
flat- and rough-bottom configurations (Figs. 9a,d and 9b,e). This is because the bottom horizontal 
flow, u, approaches zero in the strong-drag regime, and the impact of topography on QG flows 
is proportional to u- Vh, where h is the bottom topography. Our QG basin simulation results are 
consistent with the finding from previous studies (e.g., Nadeau et al., 2013) that use of realistic 
rough topography increases baroclinicity (e.g., compare upper and lower layer kinetic energies in 
their Table 2). 

It seems clear that rough topography acts to reduce the sensitivity of eddy horizontal length 
scales to bottom drag strength. Other differences between our HYCOM simulations and many 
classic studies of QG turbulence include vertical resolution (e.g., the number of layers, which 
is often only two in QG turbulence models); horizontal inhomogeneities; and other modeling 
choices, such as the choice of linear versus quadratic parameterizations of bottom drag. Although 
it is difficult to make a direct comparison, the use of a quadratic bottom drag instead of a linear drag 
may also account for part of the weakened sensitivity in HYCOM. Arbic and Scott (2008) showed 
that the sensitivities in QG turbulence to linear drag are greater than those for quadratic drag, as 
can be seen in Fig. 8 here. It seems unlikely that the reduced sensitivity seen in our HYCOM 
simulations (relative to classic studies) is strongly affected by vertical resolution, ageostrophic 


dynamics, or horizontal inhomogeneity. In support of this statement, we note that Hurlburt et al. 
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(2008) used a realistic OGCM similar to HYCOM, but with a flat bottom. They find much larger 
changes in mean SSH in response to changes in bottom drag than we see, despite the inclusion of 
horizontal inhomogeneity, ageostrophic dynamics, and higher vertical resolution in their model. 
A working hypothesis for why rough topography acts to reduce the sensitivity of eddy horizon- 
tal length scales to bottom drag is that barotropization of baroclinic energy gets short-circuited in 
the presence of rough bottom topography. Barotropization of baroclinic energy extracts baroclinic 
energy from scales near the deformation radius and injects it into the barotropic mode, typically 
at somewhat larger horizontal scales. This energy remains resident in the barotropic mode, es- 
sentially until it is removed by bottom friction. With rough topography, much of this barotropic 
energy can be transferred back to the baroclinic mode; that is, interaction between topography 
and the barotropic streamfunction forces the baroclinic mode. Assuming this to occur at a com- 
parable or faster rate than the rate at which bottom drag acts to remove barotropic energy, the 
barotropization process becomes effectively “short-circuited”. Our hypothesis and those posed 
by previous studies (e.g., Hurlburt et al., 2008) on the influence of rough topography on eddying 
flows would explain the relatively small changes observed in geostrophic SKE, SSH variance, and 


eddy horizontal length scales here. 


d. Effect of wave drag 


The strong and weak values of bottom drag used here help to demonstrate the impact of bottom 
drag strength on eddy statistics, but these extreme drag values lie outside of physically plausible 
limits. Aside from the “mid” value of Cy = 2.5 x 10-3, an additional plausible momentum sink in 
the ocean is that associated with wave drag, as described in Trossman et al. (2013, 2016). Here, we 


briefly investigate whether the sensitivity of eddy statistics to the presence of a physically plausible 
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wave drag momentum sink is qualitatively similar to the sensitivity seen with the extremes of 
bottom drag strength discussed in previous sections. 

Including wave drag and boosting bottom drag strength impact HYCOM in a qualitatively sim- 
ilar manner. The near-bottom flows are also weakened in the HYCOM simulation with wave drag 
such that the vertical profile of KE is more baroclinic relative to the simulation without wave drag 
(Fig. 10; Table 3; Trossman et al., 2016 - their Figs. 1la-d). As with the sensitivity of HYCOM 
eddy length scales to bottom drag strength (Fig. 8d; Table 5), Lxg in HYCOM is fairly insensi- 
tive to the presence of wave drag (Table 5). Area-averaged SSH variance and geostrophic SKE 
in HYCOM are both sensitive at the ~20% level to the inclusion of wave drag (Trossman et al., 
2016 - their Figs. 5 and 7; their Table 2; also see Table 2 in this paper). Lastly, the conversion rate 
between kinetic and potential energy must change with the same sign when wave drag is included 
as when bottom drag strength is increased. 

The responses of the HYCOM simulations with wave drag and strong bottom drag, however, are 
not identical. When wave drag is included, the SSH variance and geostrophic SKE are actually 
decreased, in contrast to the slight increases seen with increasing bottom drag (Table 2). This 
demonstrates the fundamentally different physical consequences of including wave drag relative 
to boosting bottom drag. Here we surmise that the spatially varying coefficient, rg,a¢, in the wave 
drag parameterization is the source of the qualitatively different responses of SSH variance and 
geostrophic SKE to the presence of wave drag as opposed to increasing bottom drag strength. From 
the results of Hurlburt and Hogan (2008), who varied bottom drag strength using only six layers 
in the vertical direction and a flat bottom in a model very similar to HYCOM, we suggest that 
applying a bottom drag over a much larger bottom layer thickness than in our HYCOM simulations 
would not cause qualitatively different behavior in the geostrophic SKE and SSH variance. We 


also suggest, based upon the horizontally homogeneous QG turbulence results of Arbic and Scott 
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(2008), that using a linear, as opposed to quadratic, bottom drag near the seafloor is not the cause 


of the qualitatively different behaviors seen when using wave drag versus bottom drag. 


5. Conclusions 


The present study investigates the sensitivity of midocean eddy statistics to bottom drag, rough 
topography, and wave drag in models of varying complexity. A primary focus is on whether 
the conclusions drawn from horizontally homogeneous, two-layer, flat-bottom, f-plane, doubly 
periodic QG turbulence simulations about sensitivity to bottom drag (e.g., Arbic and Flierl, 2004; 
Arbic and Scott, 2008) qualitatively apply to more realistic ocean models. In the QG basin and 
realistic OGCM simulations with strong bottom drag studied here, the KE is reduced in the bottom- 
most layer and generally becomes more baroclinic, as in the earlier two-layer doubly periodic QG 
results. As a result, the agreement with the vertical structure, or baroclinicity, of eddy KE in 
current meter observations is better for the OGCM simulations with a nominal “mid” value of 
bottom drag than for the OGCM simulations with a strong bottom drag. In the QG basin and 
OGCM simulations with weak bottom drag studied here, the KE becomes more barotropic, again 
in accordance with earlier two-layer doubly periodic QG results. However, the sensitivity of the 
baroclinicities in the QG basin simulations to bottom drag is reduced for rough bottom conditions 
relative to flat bottom conditions, suggesting that rough topography mediates the sensitivity of 
baroclinicity to bottom drag. 

The qualitative results about the horizontal eddy length scales seen in horizontally homogeneous, 
two-layer, flat-bottom, f-plane, doubly periodic QG turbulence damped by very weak or strong 
bottom drag are not seen in the QG basin simulations performed here with rough topography. 
In line with earlier results (e.g., Treguier and Hua, 1988), the use of rough topography reduces 


the sensitivity of eddy horizontal length scales to bottom drag strength in QG basin simulations. 


26 


Our QG basin simulations suggest that the bathymetry of the more realistic OGCM simulations 
is partially responsible for the relatively weak impact of bottom drag or wave drag on horizontal 


eddy length scales. 
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Fig. 7. 


Fig. 8. 


The rough bottom topography used in the QG basin simulations. The colorbar values are 


given in units of meters below the sea surface. The minimum depth is greater than 10 meters. 


Shown are the time-averaged geostrophic surface kinetic energies (SKE; units in m? s~*) 
in the North Atlantic, computed using a nine-point stencil (Arbic et al. 2012) from the final 
year of (a) the mid-bottom drag 1/25° global HYCOM simulation, (c) the strong-bottom 
drag 1/25° global HYCOM simulation, (b) the mid-bottom drag 1/12° Atlantic HYCOM 
simulation, (d) the weak-bottom drag 1/12° Atlantic HYCOM simulation, and (e) over all 
years (1992 — 2008) of AVISO data. 


Shown are the sea surface height (SSH) variances (units in m7’) in the North Atlantic, aver- 
aged over the final year of (a) the mid-bottom drag 1 /25° global HYCOM simulation, (c) the 
strong-bottom drag 1/25° global HYCOM simulation, (b) the mid-bottom drag 1/12° At- 
lantic HYCOM simulation, (d) the weak-bottom drag 1/12° Atlantic HYCOM simulation, 
and (e) over all years (1992 — 2008) of AVISO data. 


(a) The horizontal locations (magenta circles) of the current meter observations used in this 
study. (b) The geometric averages (solid curves) of the kinetic energy profiles over all of 
these horizontal locations. Panel b employs daily-averaged output of the strong-bottom drag 
1/25° global HYCOM simulation (red), mid-bottom drag 1/25° global HYCOM simulation 
(blue), and low-pass filtered current meter observations (black). 


Shown are the base-10 logarithms of the ratios of the kinetic energy (KE) averaged over 
the top 100 meters to that averaged over the bottom 500 meters, each computed as a time 
average over the final year of (a) the mid-bottom drag 1/25° global HYCOM simulation and 
(b) the strong-bottom drag 1/25° global HYCOM simulation. 


Shown are the base-10 logarithms of the ratios of the kinetic energy (KE) averaged over 
the top 100 meters to that averaged over the bottom 500 meters, each computed as a time 
average over (a) the final year of the mid-bottom drag 1/12° Atlantic HYCOM simulation 
and (b) the weak-bottom drag 1/12° Atlantic HYCOM simulation. 


Shown is the base-10 logarithm of the barotropic kinetic energy, KEgr (units in m? s~”), av- 
eraged over the final year of (a) the mid-bottom drag 1/25° global HYCOM simulation, (c) 
the strong-bottom drag 1/25° global HYCOM simulation, (b) the mid-bottom drag 1/12° 
Atlantic HYCOM simulation, and (d) the weak-bottom drag 1/12° Atlantic HYCOM simu- 
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Shown are results from the horizontally homogeneous, two-layer, flat-bottom, f-plane, doubly periodic QG 
turbulence simulations with linear and quadratic bottom drags from Arbic and Flierl (2004) and Arbic and 
Scott (2008); the QG $-plane basin simulations with a flat bottom and rough bottom topography; and the 
1/12° Atlantic and 1/25° global HYCOM simulations. We show non-dimensional eddy statistics: (a) the 
ratio of the domain-averaged kinetic energy (KE) in the top layer (subscript 1) to that in the bottom layer 
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ratio of the eddy length scales associated with KE in the barotropic mode (Lgr) to the Rossby radius of 
deformation (Lz), and (d) the domain-averaged ratio of the eddy length scales associated with KE in the 
upper layer (L;) to Ly. A domain average has been taken over a region (between 59.3° — 39.3°W and 
19.6° — 39.6°N) very close to the one shown in Fig. 1 for the 1/12° Atlantic and 1/25° global HYCOM 
simulations. Over this domain, Ly is assumed to be 30 km for not only the QG simulations, but also for the 
HYCOM simulations. The abscissa in each panel shows the nondimensional friction, as defined by Arbic 
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Fig. 9. 


Fig. 10. 


Shown are representative snapshots of the streamfunction (units in m* s~!) in the top layer 
of the QG f-plane basin simulations with (a-c) a flat bottom and (d-f) rough bottom topog- 
raphy. The simulations use a linear bottom drag coefficient of (a,d) 8 x 10710 5-1 (b,e) 
8x 10-8 s—!, and (c,f) 8 x 10~-© s-!. The axes have the same latitude and longitude labels 
as in Fig. 1. : 


Shown are the base-10 logarithms of the ratios of the geostrophic surface kinetic energy 
(KE) to the KE averaged over the bottom 500 meters, each computed as a time average over 
the final year of (a) the mid-bottom drag 1/25° global HYCOM simulation without wave 
drag and (b) the 1/25° global HYCOM simulation with wave drag. ye at 4 
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Fic. 1. The rough bottom topography used in the QG basin simulations. The colorbar values are given in 


units of meters below the sea surface. The minimum depth is greater than 10 meters. 
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Fic. 2. Shown are the time-averaged geostrophic surface kinetic energies (SKE; units in m? s~*) in the North 


fa) 


Atlantic, computed using a nine-point stencil (Arbic et al. 2012) from the final year of (a) the mid-bottom drag 
1/25° global HYCOM simulation, (c) the strong-bottom drag 1/25° global HYCOM simulation, (b) the mid- 
bottom drag 1/12° Atlantic HYCOM simulation, (d) the weak-bottom drag 1/12° Atlantic HYCOM simulation, 
and (e) over all years (1992 — 2008) of AVISO data. 
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890 Fic. 3. Shown are the sea surface height (SSH) variances (units in m7) in the North Atlantic, averaged over 
eo: the final year of (a) the mid-bottom drag 1/25° global HYCOM simulation, (c) the strong-bottom drag 1/25° 
see global HYCOM simulation, (b) the mid-bottom drag 1/12° Atlantic HYCOM simulation, (d) the weak-bottom 
es drag 1/12° Atlantic HYCOM simulation, and (e) over all years (1992 — 2008) of AVISO data. 
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FIG. 4. (a) The horizontal locations (magenta circles) of the current meter observations used in this study. (b) 
The geometric averages (solid curves) of the kinetic energy profiles over all of these horizontal locations. Panel 
b employs daily-averaged output of the strong-bottom drag 1/25° global HYCOM simulation (red), mid-bottom 
drag 1/25° global HYCOM simulation (blue), and low-pass filtered current meter observations (black). 
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398 Fic. 5. Shown are the base-10 logarithms of the ratios of the kinetic energy (KE) averaged over the top 100 
esa meters to that averaged over the bottom 500 meters, each computed as a time average over the final year of (a) 
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FIG. 6. Shown are the base-10 logarithms of the ratios of the kinetic energy (KE) averaged over the top 100 
meters to that averaged over the bottom 500 meters, each computed as a time average over (a) the final year of 


the mid-bottom drag 1/12° Atlantic HYCOM simulation and (b) the weak-bottom drag 1/12° Atlantic HYCOM 
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© 


0s 1/25° global HYCOM simulation, (b) the mid-bottom drag 1/12° Atlantic HYCOM simulation, and (d) the 
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FIG. 8. Shown are results from the horizontally homogeneous, two-layer, flat-bottom, f-plane, doubly periodic QG turbulence 
simulations with linear and quadratic bottom drags from Arbic and Flierl (2004) and Arbic and Scott (2008); the QG B-plane basin 
simulations with a flat bottom and rough bottom topography; and the 1/12° Atlantic and 1/25° global HYCOM simulations. We 
show non-dimensional eddy statistics: (a) the ratio of the domain-averaged kinetic energy (KE) in the top layer (subscript 1) to that 
in the bottom layer (subscript 2), (b) the domain-averaged ratio of the baroclinic KE to barotropic KE, (c) the domain-averaged 
ratio of the eddy length scales associated with KE in the barotropic mode (Lg7) to the Rossby radius of deformation (Lg), and (d) 
the domain-averaged ratio of the eddy length scales associated with KE in the upper layer (L;) to Lg. A domain average has been 
taken over a region (between 59.3° — 39.3°W and 19.6° — 39.6°N) very close to the one shown in Fig. 1 for the 1/12° Atlantic and 
1/25° global HYCOM simulations. Over this domain, Lg is assumed to be 30 km for not only the QG simulations, but also for the 
HYCOM simulations. The abscissa in each panel shows the nondimensional friction, as defined by Arbic and Scott (2008) for the 
doubly periodic QG simulations, and as defined by the relative magnitude of Cy or rgg, with respect to the control simulation, for 
the HYCOM and QG basin simulations. 
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Fic. 9. Shown are representative snapshots of the streamfunction (units in m? s~') in the top layer of the QG 
B = 


linear bottom drag coefficient of (a,d) 8 x 10-! s-!, (b,e) 8 x 10-8 s~!,, and (c.f) 8 x 10-© s~!. The axes have 


plane basin simulations with (a-c) a flat bottom and (d-f) rough bottom topography. The simulations use a 


the same latitude and longitude labels as in Fig. 1. 


Se) 


a) log;9(SKE/KE,,ottom 500) (1/25° no wave drag) 


b) logi9(SKE/KE, ottom 500) (1/25° with wave drag) 
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Fic. 10. Shown are the base-10 logarithms of the ratios of the geostrophic surface kinetic energy (KE) to the 
KE averaged over the bottom 500 meters, each computed as a time average over the final year of (a) the mid- 
bottom drag 1/25° global HYCOM simulation without wave drag and (b) the 1/25° global HYCOM simulation 


with wave drag. 
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